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Introduction 

This  report  reviews  the  progress  of  work  performed  under 
contract  on  "The  Formation  and  Subsequent  Behavior  of  Aerosols" 
between  U.  S".  Army  Chemical  Corps  and  the  Institute  of  Mathe- 
matical Sciences  (formerly  Institute  for  Mathematics  and 
Mechanics),  New  York  University.   The  period  covered  is 
November  k,    1951  -  December  31,  1953. 

Except  for  the  First  Quarterly  Progress  Report  published 
on  May  27,  1952,  it  has  been  our  practice  to  submit  at  irregu- 
lar intervals,  reports  in  the  form  of  papers  discussing  finished 
pieces  of  research.   This  has  been  a  satisfactory  practice  in 
the  early  stages  of  this  research,  when  only  certain  simple 
topics  connected  with  the  general  subject  were  considered,  in 
order  to  get  some  idea  of  results  to  be  expected  and  to  develop 
appropriate  mathematical  techniques  to  handle  more  invclved 
parts . 

The  "early"  stage  is  now  practically  over.   During  the  past 
year  we  have  obtained  certain  important  results  which  allow  us 
to  treat  the  subject  with  some  generality.   With  this  experience 
we  are  now  able  to  consider  the  more  involved  parts  of  the  sub- 
ject and  have  therefore  increased  the  staff  of  the  group  work- 
ing under  the  auspices  of  this  contract.   However,  because  of 
the  increased  mathematical  difficulty  of  the  problems  now  con- 
sidered, publishable  results  will  be  obtained  more  slowly. 
Consequently,  in  addition  to  submitting  reports  in  the  form  of 
papers,  we  shall  hereafter  revert  to  the  practice  of  publishing 
Quarterly  Progress  Reports,  briefly  reviewing  whatever  work  is 
in  progress,  even  though  it  has  not  been  completed.   This,  it 
is  believed,  will  enable  other  groups  studying  Aerosols  to  keep 
in  touch  with  progress  achieved  here. 

Three  phases  can  be  distinguished  in  the  Theory  of  Aerosols. 
The  first  is  that  of  Aerosol  Formation  which  is  discussed  in 
Part  I.   Once  a  suspension  is  formed,  drops  will  decrease  in 
size  due  to  evaporation,  or  grow  because  of  condensation  of 
vapor  on  them.   This  phase  is  discussed  in  Part  II.   Finally, 
drops  will  collide  and  coalesce  due  to  their  motion  produced  by 
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turbulence  or  falling  and  also  fall  out  of  the  air.   We  are  now 
about  to  begin  a  study  of  this  phase. 

In  the  final  section  we  list  the  staff  participating  in 
this  research  and  discuss  their  attitudes. 

Numbers  in  brackets  refer  to  references.   Items  marked  by 
a  star  (#)  were  published  under  the  auspices  of  this  contract. 


Part  I.   The  Hydrodynamics  of  Aerosol  Formation 
1 ^ 

I - 1 .   Review  of  previous  research  . 

Several  mechanisms  which  account  for  the  breakup  of  bulk 
liquid  into  drops  which  may  then  form  an  aerosol,  are  more  or 
less  understood  at  present.   The  first  of  these  —  the  forma- 
tion of  a  drop  at  the  end  of  a  capillary  tube  or  eye  dropper  -- 
which  depencT.fi  upon  gravity  in  addition  to  the  surface  tension 
and  density  of  the  liquid,  and  the  tube  diameter,  operates  too 
slowly  to  be  of  interest  in  the  present  connection.   The  second 
and  most  common  mechanism  is  the  breakup  of  a  steady  jet  cr 
stream  which  issues  from  some  kind  of  hole  or  nozzle,  as  in  an 
ordinary  atomizer.   The  hydrcdynamic  theory  of   this  type  of 
breakup  is  based  on  the  fact  that  the  unbroken  steady  flow  is 
unstable.   This  means  that  slight  disturbances  of  the  surface 
introduced  at  some  point  of  the  jet  will  grow  in  amplitude  in 
the  direction  of  flow  until  the  amplitude  of  the  surface  dis- 
turbance is  large  enough  for  opposite  surfaces  of  the  jet  to 
touch  each  other,  thus  pinching  off  the  jet  and  forming  drops. 

A  theory  of  this  type  of  breakup  was  given  by  Plateau  *■   , 

Prom  an  unpublished  review  on  "The  Hydrodynamics  of  Aerosol 
Formation"  by  J.  B.  Keller. 
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and  greatly  Improved  by  Rayleigh  [2].   The  theory  shows  that, 
for  a  jet  of  circular  cross-section,  there  is  a  certain  range  of 
wavelengthB  of  disturbance  for  which  the  disturbance  amplitude 
grows.   This  is  called  the  unstable  range,  and  in  this  range  a 
wavelength  equal  to  U.51  times  the  jet  diameter  corresponds  to 
the  greatest  growth  rate.   This  is  the  wavelength  of  maximum 
instability,  and  it  is  this  wavelength  which  accounts  for  the 
breakup  and  thus  determines  the  resulting  drop  size.   The  ampli- 

tude  grows  according  to  the  factor  e   -^         .   Here  T  is 
the  surface  tension,  p  the  density,  a  the  radius  and  v  the  velo- 
city of  the  jet,  while  z  is  distance  along  the  jet.   This  formula 
accounts  for  the  observations  of  Savart  [3]  on  the  length  of  the 
continuous,  unbroken,  portion  of  a  jet.   It  is  also  possible  to 
apply  this  theory  to  determine  the  periodic  changes  in  shape  of 
a  slightly  non-circular  jet,  originally  studied  by  Bidone,  which 
were  accounted  for  in  a  somewhat  cruder  manner  by  Magnus  [I4.]  and 
Buff  [51  •   The  theory  of  circular  jet  breakup  has  also  been 
extended  to  jets  of  viscous  liquid  and  to  jets  of  air  in  liquid 
hy  Rayleigh  [6], [7],  and  to  charged  jets  by  Basset  [8]. 

Rayleigh  [9]  has  also  given  a  two-dimensional  theory  of  the 
instability  of  a  plane  interface  between  two  fluids  in  relative 
motion,  the  fluids  being  either  bounded  or  extending  to  infinity 
in  the  direction  perpendicular  to  the  interface.   He  also  con- 
sidered the  instability  of  a  two-dimensional  jet,  and  examined 
the  effect  of  varying  the  profile  of  the  steady  flow. 

A  third  mechanism  of  drop  formation  acts  when  an  explosive 
surrounded  by  a  liquid  layer,  is  detonated.   As  Taylor  [10]has 
shown,  liquid  surfaces  may  be  unstable  when  accelerated  normal 
to  themselves,  and  it  is  this  instability  which  operates  in 
explosives  or  some  other  transient  devices.   Taylor's  theory  of 
this  instability  has  been  experimentally  verified  by  Lewis  t 11 ] - 

After  drops  have  been  formed,  they  may  possess  enough 
energy  of  vibration  or  deformation  to  break  into  still  smaller- 
drops.   A  start  on  the  theory  of  this  phenomenon  was  made  by 
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Rayleigh  [12],  and  also  by  Webb  [13],  who  determined  the  small 
vibrations  of  spherical  drops.   It  was  further  pursued  by 
Payleigh  [11+]  who  took  account  of  the  presence  of  charge  on  the 
drop  and  showed  that  enough  charge  made  the  drop  unstable. 

If  drops  are  subjected  to  a  steady  air  flow  or  to  a  shock 
wave,  they  may  also  break  up.   In  England  the  Porton  experiments 
demonstrate  these  effects,  and  rough  theories  have  been  given 
by  Penney  and  by  Taylor  in  Porton  reports. 

1-2   Breakup  of  steady  jets . 

On  the  basis  of  Rayleigh' s  work,  it  seems  that  the  breakup 
of  most  steady  streams  from  nozzles  is  due  to  instability 
governed  by  surface  tension  and  inertia.   It  should  be  possible 
to  analyze  the  breakup  of  various  fluid  sheets,  such  as  the 
"tulip"  formed  by  a  conical  nozzle,  by  appropriate  modification 
of  Rayleigh' s  theory.   In  order  to  do  this,  it  is  necessary  to 
obtain  first  an  exact  or  aporoximate  solution  of  the  steady, 
undisturbed  flow.   Since  the  exact  solution  of  jet  problems  is 
difficult,  various  approximate  solutions,  such  as  that  given  by 
the  hydraulic  theory,  will  have  to  be  used. 

In  line  with  this  program  Keller  and  Weitz  [15*-]  developed 
a  systematic  method  for  handling  thin  heavy  jets,  both  steady  and 
unsteady.  They  obtain  formal  expansions  in  powers  of  a  thickness 
parameter,  the  zero  order  terms  being  identical  with  the  hydraulic 
approximation.   The  method  is  similar  to  that  used  by  Priedrichs 
[16],  and  Keller  [17]  in  the  shallow  water  theory,  and  although 
it  is  not  known  whether  these  expansions  converge,  the  recent 

proof  of  existence  of  a  solitary  wave  in  channels,  based  on  a 

2 

similar  expansion  ,  indicate  that  they  may  converge. 

The  next  step  would  be  to  discuss  the  time  dependent  pertur- 

2 

By  Priedrichs,  K.  0.  and  D.  Hyers;  to  appear  in  the  Communi- 
cations on  Pure  and  Applied  Mathematics. 
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bations  on  the  zero  order  solutions  of  [l5-;;-]  in  a  manner  similar 
to  that  discussed  in  1-3.   1  ;e  plan  to  carry  this  out  for  the  flow 
issuing  from  a  vertical  nozzle,  and  for  the  "tulip"  flow  formed 
by  a  conical  nozzle.   This  research  will  be  continued  after  the 
completion  of  the  problems  outlined  in  1-3  now  under  consideration. 
See  also  the  discussion  in  I-lj.. 

I -3   Study  of  breakup  of  thin  liquid  layers . 

The  ideas  of  Taylor  [10]  have  been  applied  by  Keller  and 
Kolodner  [  18-::-]  in  studying  drop  formation.   They  consider  a  plane 
liquid  layer,  accelerated  in  a  direction  perpendicular  to  the 
layer  surface,  by  imposing  a  pressure  difference  on  opposite  sur- 
faces.  (The  origin  of  this  pressure  difference  need  not  be  speci- 
fied)  The  zero  order  solution  is  a  parallel  flow,  with  velocity 
of  any  of  the  two  bounding  surfaces,  which  is  also  an  exact  solu- 
tion.  It  is  then  argued  that  the  flow  deviates  from  the  zero 
order  solution  because  either  the  bounding  surfaces  are  not  per- 
fectly plane  or  the  pressure  on  the  boundary  is  not  exactly  con- 
stant.  A  first  order  perturbation  is  then  considered  and  it  is 
shown  that  to  this  order  the  surface  may  be  imagined  to  be  divided 
into  rectangles  by  two  sets  of  waves  in  perpendicular  directions. 
Alternate  rectangles  are  raised  and  lowered,  and  for  certain  dimen- 
sions the  amplitude  grows  until  opnosite  faces  of  the  liquid  layer- 
touch,  leading  to  breakup  into  drops  the  size  of  which  depends 
upon  rectangle  size.  A  certain  set  of  rectangles  leads  to  maximum 
instability,  and  determines  the  drop  size,  just  as  in  the  theory 
of  jets.   The  formula  obtained  for  the  drop  radius  r  is 

(  97iTh2    )  1//3 

\2Ip2-p-lI; 


r  = 


where , 


h  =  thickness  of  the  layer 
T  =  coefficient  of  surface  tension 
Pp-p,  =  difference  of  pressures  executed  on  opposite 
sides  . 
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The  time  of  breakdown  is  also  estimated. 

As  is  seen  from  the  formula,  the  assumption  that  there  is  sur- 
face tension  is  essential  in  this  theory.   The  theory  was  con- 
structed under  the  assumption  that  the  liquid  is  incompressible 
and  inviscid.   We  also  attempted  to  consider  the  effects  of  vis- 
cosity, but  this  led  to  complications  which  we  were  unable  to  re- 
solve.  This  attempt  was  consequently  shelved,  but  will  be  recon- 
sidered at  a  later  date. 

There  is  doubt  as  to  the  range  of  applicability  of  the  formula 
quoted  above,  and  a  more  refined  approach  to  problems  of  this  type 
is  now  under  investigation  (see  II-U).   Nevertheless,  the  simplicity 
Of  the  method  used  in  [ 18*J  warrants  its  application  to  the  problem 
of  accelerated  sphericalAcylindrical  layers. 

The  spherical  case  is  now  near  completion,  and  barring  unfor- 
seen  difficulties,  a  report  on  this  piece  of  research  will  be  pub- 
lished within  a  month  .  A  hint  on  the  method  of  solution  was  sup- 
plied in  [19*],  section  III^.  The  problem  possesses  a  great  deal 
of  similarity  to  that  of  underwater  explosion  bubbles,  and  some  of 
the  techniques  developed  for  treatment  of  the  latter  by  Keller  and 
Kolodner  [20]  are  used.  The  study  of  cylindrical  layer  will  fol- 
low that  of  spherical  layer. 

T-I4.  Rigorous  treatment  of  motion  of  liquid  layer 
We  suggested  that  the  formula  of  Part  1-3,  inasmuch  as  it  is 
of  considerable  practical  interest,  be  tested  by  exploding  into 
liquid  sheets,  and  thus  accelerating  them  into  air.   To  obtain  a 
variety  of  drop  sizes,  we  suggested  that  the  lower  surface  of  the 
layer  be  initially  deformed  in  a  controlled  way,  for  example,  by 
introducing  slight  dimples  distributed  periodically. 


3    by  Keller,  J.  B.,  I.  I.  Kolodner  and  M.  Milgram. 

^    V'ith  some  errors.   At  that  time  we  were  unaware  of  the  impor- 
tance of  surface  tension  in  problems  of  this  type. 
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Such  experiments  were  carried  out  by  Holland-  at  Camp  Dietrich, 

Maryland,  and  in  the  cases  tested  the  initial  dimples  deepened  at 
first,  whereas  at  a  later  time  small  jets  would  emanate  from  the 
troughs  of  the  dimples  in  the  direction  opposite  to  that  of  motion 
of  the  layer.   See  Pig.  la,  b,  c. 
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Oval  communication  by  Dr.  T.  Holland. 
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In  our  opinion,  these  results  do  not  necessarily  reflect  on 
the  validity  of  the  formula  of  section  1-3,  nor  do  they  indicate 
that  the  mechanisms  proposed  by  Rayleigh  and  Taylor  to  account 
for  the  breakup  of  bulk  liquid  cannot  lead  to  results  in  agreement 
with  experiment.   Our  treatment  applies  to  thin  jets  and  layers, 
but  the  mathematical  techniques  used  are  too  crude  to  account  for 
the  finer  features  observed  by  Holland,  or  to  predict  the  maxi- 
mum thickness  of  layer  for  which  our  theory  applies.   It  has  been 
suggested  that  compressibility  effects  are  responsible  for  the 
appearance  of  jets  in  the  case  tested.   This,  in  our  opinion,  is 
not  necessarily  so,  for  a  jet  will  form  whenever  two  streams  im- 
pinge on  each  other  and  this  seems  to  be  the  situation  in  Fig.  lc . 

In  any  case,  since  our  theory  does  not  explain  Holland's 
experiments,  a  better  theory  is  called  for.   The  discrepancies  are 
due  either  to  wrong  physical  assumptions  about  the  liquid  -  such 
as  incompressibility ,  or  to  too  crude  mathematical  treatment  - 
extrapolating  the  method  of  perturbations  to  account  for  large 
scale  phenomena  used  in  [l8-«-].   Accounting  for  compressibility 
effects  will  result  in  mathematical  complications  which  cannot  be 
mastered  at  present,  and  ,  in  view  of  the  argument  of  the  previous 
paragraph,  should  not  lead  to  an  essential  improvement  of  the 
theory.   On  the  other  hand,  a  vigorous  treatment  of  motion  of 
accelerated  layers  should  lead  to  a  clear-cut  answer  to  whether 
the  incompressibility  assumption  is  satisfactory,  and  if  so, 
should  shed  more  light  on  how  the  breakup  occurs. 

This  problem  is  very  complicated  and  for  simplicity  we  con- 
sider as  a  test  problem  a  two  dimensional  liquid  layer  extending 
to  infinity,  with  a  single  diaple  on  the  lower  surface  accelerated 

upwards j  see  Pig.  2. 
...  x       Placing  the  y  direction  in 
the  direction  of  accelera- 
tion, the  equation  of  the 


y-=_h+7  (x)    upper  surface  is 
Jo 


acceleration 

Figure  2 
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1    4-2 

y  =  -  at 
2 

and  that  of  the  lower  surface  is  initially 

J   =   -   h  +   zQ(x) 

We  assume  that  z   is  an  even  function  of  x,  and  that  for  x  >  0  it 
decreases  rhdnotonically  from  some  positive  value  cT  =  z  (0),  cT<  h, 
to  zero.   The  problem  is  to  find  the  boundary  function  z(x,t)(in 
which  we  are  primarily  interested),  and  the  velocity  potential 
<J)(x,y,t),  which  satisfy  the  following  equations: 

(1)  Z^4  =  0  ,  for  -at2  -  h  +  z(x,t)  <  y  <  -  at2 

2  2 

(2)  <L  =  at   for  y  =  i  at2 

y  2 

(3)  4y  -  zx4x  =  at  +  zt         I 

\  for  y  =  -^at  -h+z(x,t) 
(14.)  p(<L  +  |(V  4)2)+  p0  -  2TH  =  f(t)/ 


(5)  Iv4l  — >  0 


as   x   — >  co 


(6)   z(x,0)  =  zQ(x) 


zt(x,0)  =   0 


Here  p  is  the  density  of  liquid,  p   the  pressure  of  adjacent  air, 
T  the  coefficient  of  tension,  and  H  the  mean  curvature  of  the  upper- 
surface; 

z 

H  = XX 


2(l+z2)3/2     ' 

The  function  f(t)  is  an  arbitrary  function  of  time.   Introducing 
new  variables  y  and  <£(x,y,t)  defined  by 


-    10   - 

12- 
y  =  -   at^   +  y 

2 

4   =    aty   +  ^(x,y,t) 

the   first   four  equations   become   on  dropping  bars, 
(1')       A4>   =    0        for        -   h   +    z(x,t)    <   y   <   0 
(2«)         <j>      =   0        for        y  =   0 

OM      4y-zx4x=zt 

,  p  z  (    for  y  =  -h+z(x,t) 

(i+«  )       <L   +  -(  74)      +  az  +  i  ^£-172  =  °    J 

2  p   (l+zx)J/ 

whereas  the  last  two  remain  unchanged.   Here  we  set 

f(t)  =  p.  -  -  g2t2  +  ph   , 
°   2 

for  convenience. 

The  problem  thus  formulated  is  a  free  boundary  value  problem 
for  the  potential  equation  to  which  the  recently  developed  method 
for  handling  such  problems  may  be  applied.  This  method  consists 
in  reduction  of  the  problem  to  a  system  of  two  integro-differential 
equations  for  the  unknown  boundary  function  z(x,t)  and  for  the 
potential  o_n  the  boundary  \j/(x,t)  =  4(x,-h+z(  x,  t ) ,  t )  .  In  our  case 
these  equations  are: 


6 

By  Kolodner,  I.  I.,  to  appear  in  Communications  on  Pure  and 
Applied  Mathematics. 


-       •: 


•  :■.- 


j 


j.   . . .,   ,  y 


\  ■■' » 


• 


-  11  - 


*x-2VxZt"Zt    _ 


z 


XX 


fc        2(i+z;)  P    u+»xr' 

oo 


(7)<  *(x,t)°  i  (  ♦<*,t)iL(tan1  'lii*  >-*<*>*!+  ta-l  z(r,t)+z(x,t)-2h)d£ 
k  +  i  J  zt(r»t) {l0S[(f-x)2+(z(£,t)-Z(x,t))2][(f-x)2+(z(f,t)+ 

-  CO 

+   z(x,t)-2h)2]j  dC 

subject  to  the  initial  conditions  (6)  . 

These  equations  are  now  under  study,  but  because  of  their 
complexity  this  may  take  a  considerable  amount  of  time,  he   expect 
to  devise  on  the  basis  of  this  study  an  approximation  procedure 
which  is  more  refined  than  that  used  in  [18*-],  and  which  will 
bring  out  the  influence  of  the  thickness  of  the  layer  on  the  type 
of  behaviour  that  may  occur.   If  this  leads  to  a  theoretical  pre- 
diction in  agreement  with  experimental  results,  we  shall  have  here 
a  basis  for  predicting  the  apper ranee  of  small  jets  in  addition  to 
drops.   These  in  turn,  could  be  handled  by  the  Bay lei gh  theory 
(see  section  1-2),  and  this  would  be  a  theory  of  formation   of 
satellite  drops.   Since  the  three-dimensional  case  of  a  layer  with 
several  circular  dimples  may  be  handled  in  a  similar  manner  (al- 
though the  corresponding  equations  are  more  complicated),  it  is 
hoped  that  the  experience  thus  gained  in  considering  this  test  pro- 
blem will  be  applicable  to  the  theoretical  study  of  experiments  by 
Holland . 
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Part  II.  Decay  of  drops  by  evaporation  and  growth  by  condensation 
I I - 1 .  Review  of  pr evious  research;  formulation  of  problems 

Once  the  breakup  of  sheets  or  sprays  is  complete,  drops  will 
still  decrease  in  size  due  to  evaporation,  or  grow  because  of 
condensation  of  vapor  on  them.   Thecondensation  phenomenon  may 
also  be  used  to  produce  drops,  which  first  form  on  dust,  ions  or 
other  nuclei  and  then  grow  according  to  the  theory  to  be  dis- 
cussed below.   This  is  the  principle  of  the  Chemical  Corps'  oil 
smoke  generator,  which  is  also  used  in  artificial  production  of 
rain. 

Experimental  investigations  and  approximate  theoretical 
treatments  of  these  effects  have  been  described  by  a  number  of 
writers.   For  example,  Langmuir  [21]  and  Fuchs  [22]  treated  the 
theory  of  the  isolated  stationary  droplet,  while  Topley  and 
V'hytlaw-Gray  [23]  considered  this  problem  experimentally.   The 
evaporation  of  a  moving  droplet  was  studied  experimentally  by 
Takahasi  [21+]  and  Fr-Hs sling  [25],  and  theoretical  and  experi- 
mental results  were  given  by  Kinzer  and  Gunn  [26].   The  problem 
of  the  growth  of  a  collection  of  aerosols  or  hydrosols  has  been 
treated  extensively  by  La  Wer  and  his  associates  (see  e*g.  [27] 
and  [28]) and  by  Collins  and  his  associates  (see  e.g.  [29]  and 
[30]). 

However,  most  of  the  theoretical  work  on  this  subject  has 
been  carried  out  by  physical  chemists  who  were  forced  to  make  more 
or  less  restrictive  physical  assumptions  in  order  to  reduce  the 
problem  to  one  which  could  be  handled  by  standard  mathematical 
techniques.   For  example,  the  mathematical  difficulty  occasioned 
by  the  moving  boundary  at  the  surface  of  the  droplet  has  invari- 
ably been  avoided  by  the  introduction  of  additional  physical 
restrictions.   One  of  the  main  objectives  of  our  work  has  been 
the  exact  mathematical  treatment  of  the  moving  boundary  problem, 
as  it  stands. 


-  ■ 


> 


■  4.  - 


• 


• 


, 


1 


-  13  - 

Our  considerations  are  based  on  the  following  physical  model. 
If  a  liquid  drop  is  surrounded  by  a  gas  in  which  the  pressure  of 
the  drop's  vapor  is  below  the  equilibrium  vapor  pressure,  the 
vapor  diffuses  away  from  the  drop  surface  and  the  drop  evaporates. 
This  evaporation,  which  maintains  the  vapor  pressure  near  the  sur- 
face at  the  equilibrium  value,  results  in  decay  of  the  drop.  Con- 
versely, if  the  vapor  pressure  in  the  surrounding  gas  exceeds  the 
equilibrium  value,  vapor  diffuses  toward  the  drop,  condenses  on 
its  surface  and  the  drop  grows.   (The  same  description  applies  to 
a  drop  dissolving  or  growing  in  a  liquid,  if  the  term  "vapor 
pressure"  is  replaced  by  "concentration  of  solute".)   In  general, 
the  equilibrium  pressure  depends  on  temperature  and  on  the  radius 
of  the  drop.   In  this  model,  the  temperature  is  assumed  constant, 
and  the  dependence  of  equilibrium  pressure  on  the  drop's  radius 
is  neglected.   Hence  the  equilibrium  pressure  is  constant,  and  we 
take  it  to  be  the  equilibrium  pressure  at  a  given  temperature. 

A  more  accurate  theory  of  evaporation  (condensation)  account- 
ing for  thermal  phenomena  accompanying  this  process  is  discussed 
in  section  II-7. 

We  furthermore  assume  that  the  drop  or  drops  remain  spherical 
at  all  times  due  to  some  constraint  (like  surface  tension)  the 
nature  of  which  is  Immaterial  for  the  present  theory.   Let  a 
be  the  initial  radius  of  drop  n  centered  at  P  ,  n  =  1,  2,  3  ..., 
see  Pig.  3.   Let 

g  -  equilibrium  mass  concentration 

c  -  initial  mass  concentration  of 
o 

vapoi 
/-~Y        B  -  coefficient  of  diffusion  of 
{-<  Bn   7^  vapor  into  air 


p  -  density  of  the  condensed  vapor. 

The  mathematical  problem  is  to 
the  drop  radii   p.(t),  and  the 
mass  concentration  c(x,y,z,t)  satis- 


'-'  find  the  drop  radii   pi(t),  and  the 


Fiprure  6 
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fying  the  following  equations 


( 1 )      D  A  c   =   c . 


(2)      c(x,y,z,0)   =   c 


in  the   vapor   region 


(3)      c  =  g      on  Sn      , 


(1+)      ^(p-g)p    (t)    =   vfck  |f- 

J  J         n 


dco   , 


(5)     Pn(0)  =  an 


Here,    S      is    the    surface   of  the   n'thdroo, 
'      n 


for   all  n. 


cc 


is    the   derivative 


& 


of  concentration  in  direction  emanating  from  the  center  of  the  n'th 
drop  and  evaluated  on  S  ,  and  du)  is  the  solid  angle  element.   Con- 
dition (1+)    states  that  the  mass  of  vapor  flowing  into  the  variable 
vapor  region  is  equal  to  the  mass  of  liquid  that  evaporates. 
Introducing  dimensionless  variables  defined  by 


c(x,y,z,t)  =  cq  +  (g-c0)c(x,y,z,T) 


x  =  ax,  y  =  ay,   z  =  az 

2 

t  =  -2-t 

D 


p  ( t )  =  an  ( t )  ,   a  =  aa 

Fnv  '     Fnv  '    '    n     n 


p-g 


the  equations  become,  on  dropping  bars, 

1 


( 1 '  )       /sc 


(2»)      c(x,y,z,0)    =    0 


. 


/ 


in   the   vapor   region 


.   ,    • 


.  .    ■ 


- 
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(3')    c  -  1     on  Sn       ^ 

^)   Pn  =  fedlF-d^     L   ^  all  n. 
(5' )    P  (0)  =  a 

If  g  >  c   ,  the  drops  evaporate  and  a  >  0  .    If  g  <  c   ,  drops 
grow  and  a  <  0.   In  all  practical  cases,  g  «  p,  c   «  p,  so  that 
a  is  a  small  number  ( a  ~  2  x  10   for  water  drops  evaporating  int 
dry  air  at  73°  F) •   In  any  case,  c   <  p,  so  that  a   >  -1  . 

In  case  of  a  single  drop,  c  becomes  a  function  of  distance 
from  the  center  of  the  drop,  r.   Choosing  for  the  scaling  factor 
a"  the  initial  radius  of  the  drop,  the  equations  simplify  to 


o 


tt  „  M 


(1")  ^c  =  i  (rc)rr  =  ct   ,     r>  p(t) 

(2")  c(r,0)  =0,     r  >  1 

(3")  c(P(t),  t)  =  1   ,      t  >  0 

(V)  p(t)  =  acr(p(t),  t)   ,      t  >  0 

(5")  p(0)   =   1 

It  is  worthwhile  to  consider  also  the  case  of  one-dimensional 
"drop",  i.e.  that  of  evaporating  slab  of  water.  In  this  case,  the 
equations  become 

(l»t)  Ac  =  cxx  =  et   ,    x  >  P(t) 

(2"' )   c(x,0)  =  0  ,     x  >  1 

(3"')   c(p(t),  t)  =  1  ,   t  >  0 
(U"»)   p(t)  =  o.cx(p(t),  t)  ,   t  >  0 
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(5"')   p(o)  =  l 

with  the  additional  condition,   c(x,  t)  =  c(-x,  t)  . 

We  shall  refer  to  these  problems  as  problems  A,  B,  and  C. 

The  three  problems  formulated  above  (in  decreasing  order  of 
difficulty)  are  free  boundary  value  problems  for  the  heat  equation. 
Problems  of  similar  nature  arise  in  the  study  of  change  of  phase, 
and  were  first  considered  by  Stefan  [31],  who  solved  a  very  special 
case  explicitly.   Stefan-like  problems  were  considered  later  by 
Datsev  [32],  Evans,  Isaacson  and  KacDonald  [33],  and  Evans  [3i|],  but 
their  methods  and  techniques  turned  out  to  be  of  little  value  in  the 
study  of  problems  arising  in  cases  considered  here. 

I I - 2   Approximate  treatment  of  the  problem  of  a  single  d r o p . 

The  study  of  behavior  of  a  single  drop  is  preliminary  to  that 
of  collection  of  drops.   The  single  drop  problem  is  in  itself  of 
great  interest.   Experimental  data  concerning  evaporation  are  avail- 
able and  these  indicate  that,  except  for  a  short  initial  period  and 
very  small  drops,  the  square  of  the  drop's  radius  decreases  linearly 
witn  time;   see  [2.1+1,    [25],  [26],  Godsave  [35]  and  Ranz  and  Marshall 

[36]. 

7 
The  initial  study  was  reported  in  [19*],  section  II.   The  one- 
dimensional  case,  problem  C,  can  be  solved  explicitly  and  the  re- 
sult is 

(1)    p(t)  =  1  -  2z/F" 

where  z  is  the  solution  of  the  transcendental  equation 

~    oo     - 

-z2   ('     -a2 
z   =  2ae  e  °   do"   . 

J-z 

This  equation  has  one  solution  provided  that  o  >  -1  (physically,  the 


7 

by  Keller,  J.  B, 
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case  a  =  -1  corresponds  to  the  critical  point  when  the  change  of 

a 


phase  is  not  well  defined).   If  a  is  small,  z  f^>— —    ,  and  in  case 


of  evaporation,  a  >  0,  this  yields  for  the  time  of  evaporation  t  , 

(1.)  to  -a, 

although  equation  (1)  is  exact,  approximation  methods  were  also 
devised  with  intent  of  generalization  to  the  three-dimensional  case, 
problem  B,  where  an  explicit  solution  cannot  be  constructed.   The 
first  of  these  consists  in  assuming  an  approximation  to  the  evapora- 
tion curve,  and  foregoing  the  second  free  boundary  condition,  eq, . 
(1-li" ').   One  gets  thus  an  ordinary  boundary  value  problem  to  solve 
and  one  uses  this  solution  in  eq.  (l-V)  to  get  an  improved  approxi- 
mation for  the  evaporation  curve.   If  one  takes  for  the  first   approx- 
imation to  p 

p  =  constant  =  1 

(which  corresponds  to  a  =  0),  the  second  approximation  yields 

P  =  1  .  2o_  /F- 

/*- 

hence  a  result  identical  with  (1)  for  small  a.   Iteration  of  this 
process  leads  to  rather  complicated  expressions. 

The  second  method  is  the  so-called  quasi-static  approximation. 
One  assumes  that  the  concentration  is  slowly  time  dependent  so  that 
the  term  c      in  the  differential  equation  can  be  at  first  omitted. 
One  thus  constructs  a  first  approximation  to  the  concentration  and 
the  evaporation  curve.   The  exact  result  in  the  one-dimensional  case 
shows  however,  that  c.  is  infinite  for  t  =  0  at  x  =  1.   Thus  this 
method  is  not  expected  to  be  satisfactory  for  small  t.   Indeed  in 
this  case,  it  cannot  be  applied  at  all,  for  it  leads  to  an  expres- 
sion for  c  which  does  not  remain  bounded  as  x  increases  to  infinity. 
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However,  both  methods  can  be  applied  in  the  three-dimensional 
case,  problem  3.   The  first  yields 


for  a  small 


whereas  the  second  gives 

(3)   p2(t)  =  1  -  2ct   ,       (3')  t  ~i- 


o 


2g 


The  two  methods  thus  lead  to  different  results  and  predict  evapora- 
tion times  differing  by  a  factor  of  2  for  small  o.   Both  results 
are  therefore  inconclusive,  although  the  first  one  agrees  with  the 
exact  result  in  the  one-dimensional  case,  whereas  the  second  is  in 
agreement  with  the  experimental  results  in  three-dimensional  case. 
A  third  approximation  method  was  devised  by  T'eller,  Kolodner, 
and  Ritger  [37---].   They  represent  the  concentration  by  an  integral, 
the  form  of  which  is  suggested  by  the  exact  solution  of  the  one- 
dimensional  problem.   This  integral  satisfies  the  heat  equation  no 
matter  what  its  kernel  is,  provided  that  it  has  the  required  form. 
Kext  the  kernel  "f"  is  expanded  in  formal  powers  of  the  variable 
on  which  it  depends, 


co 

>..  n 


f(T  )  =  A0(l  +y      AnTn) 
n=l 


whereas  the  unknown  radius  is  assumed  to  have  an  expansion  of  the 
form 

p(t)  =  1  +  2  y     antn/2  , 
fel 

and  the  coefficients  are  determined  by  satisfying  the  two  free 
boundary  conditions.   As  is  shown,  the  two  expansions  can  have  at 
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best  an  asynrototic  character  neai  t=0  whereas  an  expression  for  p(t) 
valid  duiim  the  full  evaporation  period  is  desired.   To  get  this, 
one  defines  as  the   'th  Qpproxiraatior  to  f(T),  t  e  oolynomiel 
obtained  from  the  series  on  dropping  ^11  tcr.^s  frcn  the  ,+1'st  on, 
with  coefficients  ..  (n=C ,  1, . . . .  )  previously  determined,  and  one 
defines  the   ' th  aooroxiftation  to  r,  by  requiring  that  it  satisfies 
the  second  boundary  condition  (1-[|.")  with  c  replaced  by  its  h '  th 
approximation.   At  each  step  one  thus  gets  an  ordinary  differential 
equation  which  must  be  integrated  numerically.   A  calculation  was 
carried  out  for  a=2  x  10   ,  N=l,2,3,  showing  that  the  second 
approximation  for  p  agrees  with  the  third  to  3  decimals.   For  a 
small  these  equations  can  be  simplified  without  introducing  an 
appreciable  error,  and  one  gets  for  the  second  approximation, 

(h)  P2-  1  -  ^/t~  -   2at   , 

also  agreeing  to  within  3  decimals  with  the  computed  result.   The 
evaporation  time  is 


(U«)  t  ~± it i_  +  od)     . 

O  o  m —      / 


2a        3/21?  /a~ 

Formula  ([).)  agrees  with  experimental  results  over  the  full 
range  of  evaporation  time. 

I I -3   Rigorous  treatment  of  an  evaporating  drop . 
Of  the  three  methods  discussed  in  the  preceding  section,  the 
third  cannot  be  easily  adopted  to  handle  more  complicated  problems, 
whereas  the  first  two  can.   Since  the  third  method  led  to  a  result 
in  agreement  with  that  of  the  first  method  for  small  t,  and  with 
that  of  the  second  for  large  t,  its  justification  would  at  the  same 
time  establish  the  respective  ranges  of  applicability  of  the  first 
two  methods,   kith  this  in  mind,  we  attacked  the  problem  B  rigor- 
ously hoping  to  reestablish  the  formula  (2r-[j.)  . 

Since  equation  (2-1")  leads  to  a  very  simple  ordinary  differen- 


r" 


V 


. 


.   . 
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tial  equation  for  the  Laplace  Transform  of  u(r,t)  =  rc(r, t)  with  res- 
pect to  t,  use  of  the  Laplace  Transformation  seemed  to  be  proper  for 
a  problem  of  this  type.   Evans,  Isaacson  and  Tic  Donald  [33],  expressed 
the  hope  that  this  technique  should  produce  an  elegant  solution  to 
their  problem,  but  got  nowhere.   Likewise,  in  case  of  our  problem  the 
use  of  Laplace  Transform  led  to  a  formal  expression  for  the  concen- 
tration function,  which  would  be  the  correct  expression  if  the 
boundary  p  were  known,  but  could  not  be  easily  used  to  determine  their 
boundary.   (It  was  realized  much  later  that  this  expression  could  be 
obtained  in  a  simpler  way  by  applying  Green's  theorem  to  the  function 
u,  see  [38]*)   What  evolved,  however,  from  an  interpretation  of  this 
formal  result,  was  a  basically  simple,  but  very  general  and  powerful 
method  for  solving  the  free  boundary  value  problems  for  the  heat  equa- 

tion,  by  Folodner  [39]  .   The  method  consists  in  a  reduction  of  the 

q 
problem  to  an  integro-differential  equation  for  the  boundary.    This 

method  was  applied  to  the  case  of  evaporation  of  a  drop  by  Kolodner 

[1|0]    and  yields  p(t)  as  the  solution  of 


z(pfO) 


-T2 


/H   pp=  -2a  [  (1-  -l)e-T  dt+  I  £ilMHz(p,o-)e-z  <P'tf>d<j-  =  G[p] 


-00      ^ 

p(0)   =   1 
where 

8(Xfff)  =  ZzEi£L  ,    Z(x,0)  =  2£i   . 

The  reference  is  to  an  abstract  of  a  paper  presented  at  the  April 
1953  meeting  of  the  American  Mathematical  Society  in  New  York  City. 
Details  will  appear  in  Communications  on  Pure  and  Applied  Mathematics. 

9 

The  fundamental  idea  is  the  same  as  the  one  used  by  the  author  for 

free  boundary  value  problems  for  the  potential  equation,  mentioned  in 

section  I-1+,  which  actually  was  developed  later. 

10 

An  extract  from  a  longer  and  rather  involved  paper  to  appear  in 

Communications  cm  Pure  and  Applied  Mathematics. 
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This  equation  shows  that 

(2)  p2(t)  >  1  -  ^2.   yT  -  2a    , 

whereas  an  independent  argument  yields  directly  from  the  formulation 
of  the  problem,  an  upper  bound  for  p 

(3)  p2  <  1  _  2a..  t   , 

1+a 

see  Pig.  I4..   Prom  these  one  gets  a  two  sided  inequality  for  the 

evaporation  time  t   , 
e  o  ' 

(k)  — +  ^(i-/mr)<  t   < i-+  i 

bound  2a  71      ca         2a  2 

For  a=2  x  10   ,  the  case  considered 

previously,  either  of  the  bounds 

2 
(2)  and  (3)  yield  p  with  an  abso- 

y  lute  error  not  exceeding  .0072  over 

the  full  evaporation  time,  whereas 

-.      .  (1+)  determines  t   with  a  relative 

Figure  4  ^ ■  o 

error  of  at  most  .0072. 
p 
The  lower  bound  for  p   is  identical  with  formula  (2-1).).   Thus 

the  results  previously  obtained  are  confirmed,  provided  that  the 
solution  of  eq.  (1)  exists.   The  existence  of  a  unique  solution  in 
the  small  can  be  demonstrated  using  standard  techniques,  provided  that 
0  <  a  <  .57  ;   since,  however,  we  are  interested  in  results  for  the 
full  evaporation  period,  such  a  proof  is  not  satisfactory.   An 
existence  and  uniqueness  theorem  in  the  large  was  constructed  by  show- 
ing that  the  solution  of  eq.  (1)  is  the  limit  of  an"alternating" 

sequence  of  functions  p  (t)  defined  over  the  interval  0  <  t  <  t 

^n   '  —   —  n 

increasing  with  n  and  tending  to  t  as  n  — >oo.  Any  two  consecutive 
terms  of  the  sequence  form  an  upper  and  lower  bound  for  the  solution 
in  the  interval  over  which  these  terms  are  defined.   Thus  there  is  a 
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possibility  of  improving  the  bounds  (2)  and  (3)  which  is  useful 
when  either  o  is  larger  (but  smaller  than  .38),  or  when  higher  ac- 
curacy is  desired.   p   is  identical  with  the  right  hand  side  of  in- 
equality (2),  whereas   p,  is  an  upper  bound  which  is  better  than 
eq.  (3)  for  small  t,  but  not  as  good  as  (3)  for  large  t. 

1 1  -\\     Growth  of  drop  by  condensat  ion 
The  case  of  condensation  differs  from  the  case  of  evaporation 
in  that  -  1  <  a  <  0,  (see  section  II-l).   The  general  method  discussed 
in  [39]  leads  to  the  same  equation  for  p(t)  as  equation  (3^-1)  with 
a  <  0.   However,  as  is  indicated  in  [39],  the  condition  that  p(t) 
satisfies  this  equation  is  a  necessary  condition  for  the  solvability 
of  the  free  boundary  value  problem  in  case  of  condensation,  whereas  the 
sufficiency  has  not  been  proved.   Nevertheless,  it  is  possible  to 
derive  a  second  equation  for  p(t),  the  solvability  of  which  is  fully 
equivalent  to  the  solvability  of  the  problem  as  posed  originally. 
This  equation  is 


(1)    /i   p(t)    =   2p(t) 


0?  p 

'    e""(    dt-2/t   e"z    (P.°)_   1 
"  a. 

(p,0)  0 


PXcnP(o-)e-z2(p,q)dcJ> 
/xT~<s 


The  techniques  used  in  solving  equation  (3-1)  do  not  apply  in  this  case 
and,  so  far,  this  equation  has  not  been  solved  in  the  large.   Yet, 
it  can  be  shown  quite  independently,  that 

(2)  p2:.l_2o_t,-l<a<0   . 

1+a 

Probably  for  t  large,  p-W  TZ      7~ 

1+a 
We  are  now,  in  conjunction  with  the  study  of  other  problems  dis- 
cussed in  section  II-5,  in  the  process  of  devising  new  techniques, 
which,  it  is  expected,  will  enable  us  to  solve  equation  (1),  and  in 
particular  will  establish  an  upper  bound  in  the  large  for  p  . 
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I I - 5  Approximate  treatment  of  collection  of  drops 
An  approximate  treatment  of  collection  of  evaporating  drops 
based  on  the  quasi-static  approach  was  attempted  by  Keller 
Difficulties  were  encountered  which  could  not  be  overcome  without 
Introducing  certain  additional  assumptions  concerning  the  model.   As 
it  appears  now,  on  the  basis  of  a  more  accurate  treatment,  (as  yet 
incomplete),  discussed  in  section  II-6,  these  assumptions  are  un- 
necessary, but  the  method  must  be  modified.   This  work  will  be  taken 
up  again  in  the  future. 

Another  approximate  approach  for  a  colDection  of  initially 
identical  drops  with  centers  at  vertices  of  a  cubic  lattice  stems 
from  the  consideration  of  a  one-dimensional  collection.   In  this  case, 
Pig.  5,  the  evaporation  curves  are  all  the  same,  and  the  concentration 
a  i  c  is  a  periodic  function  of  x  with 

i  oeriod  2M  =  distance  between  two 

\    '     / 


\     :     M       successive  drops,  and  has  the 

/'i)\   !   /©\   !    /®\      lines  x  =  M'  x  =  3^'  •••* 

■*-  -4~     ^  >  x  =  (2n+l)M...   for  lines  of  sym- 


metry.  Clearly,  on  these  lines, 
c   =0.   To  solve  the  problem,  it 
--if^re  5  is  thus  sufficient  to  consider  a 

half  period  interval,   c  and  p  then  satisfy   equations  (1-1"')  through 
(1-5"'),  except  that  x  is  restricted  to  the  range  p(t)  <  x  <  M,  and 
for  x  =  M,  c   =  0  .   This  condition  also  states  that  the  flux  of 
vapor  through  the  line  of  symmetry  is  0,  and  therefore  the  problem  of 
a  one-dimensional  collection  is  fully  equivalent  to  that  of  a  single 
drop  placed  in  an  impermeable  box,  with  walls  at  distance  M  from  the 
center  of  the  drop. 

In  the  three-dimensional  case  the  situation  is  similar.   However, 
to  be  able  to  apply  the  general  method  of  [39],  we  assume  for  simpli- 
city that  the  box  is  not  cubic,  but  spherical.   This  assumption,  it 
is  believed,  ought  to  have  little  effect  on  the  result. 

The  method  of  [39]  leads  to  the  following  equations  for  the 
radii  : 


'    ! 


■1 


7  & 
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2,      -,      t    *,  >s  _2, 


(1)  /Vi   P(t)=   -  — e"Z    (p'0)+  (    £i£iz(p,d)e-z    (P'^dd- 


-  ( 
A  4  t-cr 


7 


2,^, 


t    . 


+  a_  e-z    (2M-P,0)_  f   £ioIz(2M-p,(5)e-z    <2M-P»*>d<r 


o 


t-o 


V.  p(o)     =     l 


in  the   one-dimensional   case,    and 

t 
(2)     (fR   P(t)p(t)    =   G[p]    +   2a  £(<J)hv(p(t),t,<J)dff 

0 
p(O)      =      1 


x' 


where,    G[p]    is   defined    in   equation    (3-1), 

t 


«(t)    =    (A-  JL)e-z    (M.O)+i.   !     p(<y)P(tfl(14lzlM££l)e-z    ^'^do" 


2/t 


2a  JQ      t   -   d  /t^cT 


h(x,t,(r): 


M/i 


:-M    |    /t^O 


exp(-( =r)     )+  -exp( +  — 2) 

/t^d  2/t-a-  M  M  M 


*/Z=3       M      e-7    dT 


-  oo 


z(p,<y)    =    z(p(t),<f), 


and  z(x,C)  is  defined  in  section  II-3.  These  equations  apply  to  the 
case  of  evaporation  only,  but  can  be  used  in  case  of  condensation  in 
the  sense  explained  in  section  II-li. 

These  equations,  as  well  as  eq.  (I4.-I),  and  some  other  equations 
(such  as  equations  (6-i|) ) stemming  from  related  problems,  are  now 
being  studied,  and  our  aim  is  to  produce,  in  each  case,  a  solution 
in  the  large  by  means  of  a  single  iteration  sequence,  which  would  also 
yield  simple  lower  and  upper  bounds  to  the  solution.   Such  bounds, 
due  to  the  small  value  of  a,  are  expected  to  be  sufficiently  sharp 
for  all  practical  purposes,  thus  eliminating  the  need  for  costly 
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numerical  integrations.  This  was,  of  course,  achieved  in  the  case 
of  a  single  droplet.  This  phase  of  our  work  is  well  advanced,  but 
the  methods  developed  for  handling  the  above  equations  are  too  in- 
volved to  be  discussed  here  with  any  detail. 

I I -6  Collection  of  drops  -  general  case 

The  method  of  [39]  has  been  recently  generalized  to  account  for 
the  free  boundary  value  problem  for  the  three-dimensional  heat  equa- 
tion, with  any  number  of  boundary  conditions.   It  is  thereby  possible 
to  reduce  the  general  problem  of  a  collection  of  drops  (problem  A) 
to  that  of  solving  a  system  of  integro-diff erential  equations;  the 
number  of  these  equations  is  equal  to  twice  the  number  of  drops. 

In  addition  to  terms  defined  in  connection  with  the  formulation 
of  problem  C,  section  II-l,  we  define; 

E,      -  position  vector  of  point  P   (Pig.  3) 

r  ,  9   ,  4   -  polar  coordinates  of  a  point  with  respect 
n   n  '  Tn 

to  a  system  of  coordinates  centered  at  r 
J  n 

with  cartesian  axes  obtained  by  translation 
of  x,  y,  z  axes. 
(9  ,  4>  )  -  unit  vector  with  origin  at  F 

do  -  solid  angle  element  on  sphere  centered  at  r 

«■  =  £(©n4nt'0  =  .5n+^nPn(f)'  P°lnt  on  sphere  Sn(f)  centered  at  I"n, 

of  radius  p  (T ) 

x        -     position  vector 


u    (t,    9    ,    <j>    )    =  -^-  c 
nv    '      n'    vn'        •>. 
or 
n 


x  =   f(9n,    cj>n,    t) 


The  equations  for  the  radii  of  the  drops  are: 


-un(t,e  ,<L)  =  -A-  r-  v 

n  m 
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*  -  l(en'^n't) 


n  =  1,2 , . . . 


(1) 


pn(t)  =f-$un(t,9n,4n)d<Jn 

,  p  ( 0)  =  a   , 

where 

<r  -a  \       /r  +a     1 

m     r_  v    %   \   2/t  2/tT     / 


Except  for  the  special  case  of  drops  of  initially  equal  radii, 
centered  at  the  vertices  of  a  cubic  lattice,  when  equations  (1)  are 
the  same  for  all  n,  we  have  here  a  large  number,  possibly  infinite, 
of  equations.   To  consider  these  equations  with  full  rigor  is  a 
hopeless  task,  and  an  exact  solution  is  not  even  desired  considering 
that  in  study  of  aerosols  Initial  data  on  individual  drops  cannot 
and  need  not  be  specified. 

A  first  simplification  of  system  (1)  is  obtained  on  the  basis 

of  a  rather  reasonable  assumption  that  the  distance  d   between  any 

nm  J 

two  drops  n  and  m  is  large  compared  to  their  radii.   It  then  turns 
out  that  after  integrating  the  first  set  of  equations  over  the  sur- 
face of  the  n*th  drop,  and  dropping  terms  which  can  be  neglected 
according  to  our  assumption,  u  (t,9  ,4>  )  can  be  eliminated  using  the 
second  set  of  equations.   We  get  thus  a  set  of  equations  of  the  form 
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(3)  J  ir 

n=l,2. . . 

TvPn(0)   =   an   . 

Here,  G[p  ]  is  a  functional  of  p   identical  with  the  right  hand  side 
of  eq.  (3-1),  whereas  the  other  terras  represent  the  corrections  to 
eq.  (3-1)  due  to  the  influence  of  other  drops. 

Equation   (3)  bears  some  similarity  to  equations  obtained  by- 
Keller  in  trying  to  solve  the  problem  of  a  collection  of  drops  by 
a  quasi-static  method,  and  at  the  same  time  explains  why  the  immodi- 
fied  quasi-static  approach  fails  to  lead  to  a  coherent  theory  of 
evaporation,  as  noted  in  section  IT-5«   As  may  be  shown,  the  order 
of  magnitude  of  f  [  pn,  Pm,an,  a^d^]  for  large  d^  is 

d2 

1     nm 

f  -  d   e~  ,  . 

nm     l+t     , 

hence,  the  series  in  (3)  converge  if  t  is  finite.   As  we  found  al- 
ready, the  quasi-static  approach  is  valid  near  evaporation  time, 
hence  for  t  large;  see  sections  II-2  and  It-3.   In  the  process  of 
carrying  out  the  simplification,  t  is  replaced  by  infinity  in  some 
places.   If  this  is  done,  without  further  scrutiny  in  the  series  in 
eq,  (3),  the  series  becomes  divergent,  for  its  terms  become  asymp- 
totic to  d    .   This  is  the  difficulty  that  was  actually  encountered 
by  Keller. 

If  the  drops  are  initially  of  the  same  radius,  and  their  centers 
are  at  the  vertices  of  a  cubic  latice,  p  =  p  =  p(t),  any  n  and  m, 
for  all  t.  Eq.  (3)  reduce  then  to  a  single  equation: 

>(t)p(t)  =  G[p]  -  G*[p]  +  V"  flp.p.l.l^^]. 

m 
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This  is  the  case  considered  already  in  section  II-5>  hut  the  ap- 
proximation is  of  a  different  nature,  and  should  produce  a  more 
accurate  result  for  p  than  eq.  (5-2).   Eq.  (k)    is  of  a  form 
similar  to  those  encountered  already  in  preceding  sections,  and 
is  now  under  investigation. 

As  to  the  general  case,  the  system  (3)  is  still  too  complicated 
for  practical  considerations.   A  further  simplification  can  be  ob- 
tained from  the  study  of  the  equation  determining  the  evaporation 
curve  of  a  single  drop,  eq.  (3-1).   On  this  basis  one  can  show  that 
certain  terms  in  (3)  can  be  so  simplified,  without  introducing  an 
appreciable  error  in  the  final  result,  that  the  system  (3)  becomes 
a  system  of  ordinary  differential  equations.  This  system  is  now 
being  studied,  and  although  we  have  no  concrete  results  so  far,  it 
appears  that  the  best  approach  would  be  that  used  in  statistical 
mechanics,  since  tut*  .system  studied  is  formally  the  same  as  the 
system  of  canonical  equations  of  motion  of  large  number  of  particles, 
with  the  p  (t)  being  interpreted  as  displacements.   If  this  line  of 
attack  proves  to  be  successful,  the  result  will  be  in  form  giving 
the  relative  number  of  drops  of  radii  p(t)  at  time  t,  in  terms  of 
the  initial  density  of  radii.   A  result  of  this  form  is  precisely 
what  is  desired  in  practical  application. 

We  have,  throughout  this  section,  spoken  exclusively  of  evapo- 
ration, but  our  considerations  apply,  with  some  modifications,  to 
the  case  of  condensation. 

I I -7.   Effects  of  heat  flow  on  process  of  evaporation 
In  the  above  treatment  of  the  evaporating  drop  (sec.  11-2   and 
II-3)  the  effects  of  temperature  have  been  neglected  in  order  to 
simplify  the  discussion  of  the  central  mathematical  difficulty 
encountered  in  finding  the  boundary  of  the  evaporating  drop,   how, 
however,  using  the  techniques  already  developed,  it  is  possible  to 
attack  the  more  general  problem  including  the  temperature  effects. 

11 
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Vo  final  results  have  been  obtained  as  yet,  but  we  shall  outline 
our  approach  and  our  plans  for  future  work  on  this  aspect  of  the 
problem. 

In  general,  there  is  a  decrease  in  temperature  in  the  vicinity 
of  the  surface  of  an  evaporating  drop  because  of  the  absorption  of 
the  latent  heat  of  vaporization.   This  results  in  a  flow  of  heat  in 
the  interior  of  the  drop  and  in  the  vapor  region  outside  the  drop. 
Moreover,  the  variation  of  temperature  produces  a  change  in  the 
equilibrium  pressure  at  the  drop  surface  and  hence  affects  the  rate 
of  evaporation. 

Several  possible  formulations  of  this  evaporation-temperature 
problem  have  been  suggested,  but  the  one  that  seems  to  lend  itself 
most  readily  to  treatment  is  one  similar  to  that  proposed  by  Kinzer 
and  Gunn  [26]. 

We  consider  an  evaporating  drop  as  in  sec.  II-l,  but  we  now 
assume  that  the  equilibrium  concentration  at  the  drop  surface  is  a 
known  function  of  the  temperature.   We  further  assume  that  the  tem- 
perature q(r,t)  is  initially  constant  everywhere,  is  continuous 
across  the  drop  boundary,  and  satisfies  the  equation  of  heat  conduc- 
tion outside  the  drop,  but  is  uniform  (a  function  of  t  only)  inside. 

The  mathematical  formulation  of  this  problem  is  as  follows  (in 
dimensionless  coordinates  introduced  by  transformations  similar  to 
those  in  II-l)  : 

(1)  rct(r,t)  =  (rc(r,t))rr  ,  rqt(r,t)  =  (rq(r,t))rr  ,  for  r>p(t),t>0. 

(2)  c(r,0)  =  0   ,   q(r,0)  =   0 

(3)  c(P(t),t)  =  F(q),  for  0  <  r  <  p(t),  t  >  0 

(k)         cr(p(t),t)  *  [1  +  |  -  P(q(p,t))]p(t),  qr(p(t),t)=  -6p(t),for  t>0, 

(5)  c(oo,t)  =  0,  q(oo,  t)  =  0,  t  >  0  , 

(6)  p(0)  =  1   ,   P(0)  =  1   , 

where  r  is  the  radial  distance  from  drop  center,  t  the  time,  q  the 
temperature,  c  the  vapor  concentration,  P  the  known  (empirical) 
function  giving  equilibrium  concentration  in  terms  of  temperature, 
p  the  drop  radius  and  a  and  6  positive  constants.   For  convenience 


> 


. 


i 


>        :    :   .    ; 
j"      :  • :  ■ 


' 


< 


' 


_  K  } 


■   ..   :        ■       :       -. 


!.  -•' ! .  U 


I     .->■  - 


:  ^';- 


-  30  - 

we  set  the  diffusion  coefficient  and  the  coefficient  of  thermal 
diffusivity  numerically  equal;  an  assumption  of  this  type  is  un- 
justified and  unnecessary  for  the  continuation  of  our  work  and 
will  be  corrected  in  the  final  report  on  this  piece  of  research. 
This  problem  is  attacked  by  the  method  of  [39].   We  obtain 
a  pair  of  simultaneous  integro-diff erential  equations  for  p't) 
and  q(p(t),  t);    namely: 

(7)     \   l/5T   p(t)P(t)    =    G[p]    +  JT  H[p,    q] 

q  =     K[p,    q   ] 


where 


H[ 


G[p]  =  right  hand  side  of  equation  (3-1) • 

i      f  r(p(<mi-F])' 


/~l    -   (f 


p,q]    =    a  -i(l-pp)(l-F(q))    +  i-      [    f 
l-  /*"      o 

+     £ e    z   do" 

t   -   C  J 

K[P,    q]    =    5pp    -  i-     \      \W*W    +    6p(cT)p(^q 

l/tT"    J       i/t~-  cr  t   -  o~ 

=      P(t)    -    p(cf) 
2/   t    -   cr 


e    z   dc~ 


These  equations  are  now  being  studied  and  we  expect  that  a 

solution  can  be  obtained  by  an  extension  of  the  iteration  pro- 
cedure which  was  used  in  solving  eq.  (3-1) . 

Assuming  that  we  can  obtain  a  solution  to  this  problem,  we 
should  then  be  able  to  take  account  of  temperature  effects   in 
the  problem  of  the  collection  of  drops  (by  means  of  the  enclosed 
drop  model,  section  II-5,  or  more  accurately,  as  in  section  II-6). 
Also,  we  may  be  able  to  extend  the  method  to  include  heat  conduc- 
tion in  the  interior  of  the  drop  or  some  of  the  other  temperature 
effects  that  have  been  neglected. 
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General  Remarks . 
1 .   Personnel 

The  following  are  now  actively  participating  in  the  research 
work  on  this  project,  either  on  part  time  or  full  time  basis: 

Dr.  E«  D.  Corson, 
Dr.  B.  Epstein, 
Dr.  I.  I.  Kolodner 

Mr.  M.  Milgram 
Mr,  P.  D.  Ritger 

It  is  planned  to  enlist  the  services  of  another  research 
assistant  within  two  months. 

The  work  is  coordinated  by  Dr.  Kolodner,  under  the  general 
supervision  of  Dr.  R.  Courant,  Director  of  the  Institute  of  Mathe- 
matical Sciences,  K.  Y.  U.   Other  members  of  the  Institute  of 
Mathematical  Sciences  are  frequently  consulted.   Dr.  J.  B.  Keller, 
on  leave  of  absence  from  New  York  University,  at  present  head  of 
the  Mathematics  Branch,  Office  of  Faval  Research,  V'ashington  D.  C., 
who  directed  the  research  on  this  contract  until  September  1953, 
continues  to  take  active  interest  in  our  work  and  serves  in  an 
advisory  capacity.   It  is  hoped  that  by  the  end  of  Summer  1951+  he 
will  rejoin  this  group  as  an  active  member. 

Other  members  of  the  Institute  who  until  recently  devoted  part 
of  their  time  to  this  project  are 

Mr,  M.  Jordan 

Ir.    G.    B.    Vi/hitham 
Dr.    B.    Zumino 
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2.  Attitudes 

It  is  realized  by  the  members  of  this  group  that  the  sponsoring 
agency  is  mainly  interested  in  results  in  the  form  of  workable  form- 
ulae and  discussions  pertaining  to  the  theoretical  description  of 
the  Formation  and  Behavior  of  Aerosols.   With  this  in  mind,  in  the 
initial  stage  of  our  work  we  attacked  various  problems  by  approxi- 
mation methods  without  paying  too  much  attention  to  mathematical 
rigor.   Correspondingly,  results  were  often  crude  and  not  easily 
justifiable. 

As  it  turns  out,  the  problems  encountered  are  of  a  complex 
mathematical  nature.   To  treat  them,  various  simplifying  physical 
hypotheses  are  made,  and  then  mathematical  approximations  follow.   If 
the  theoretical  results  thus  obtained  do  not  agree  with  experimental 
data,  it  is  hard  to  determine  whether  the  discrepancy  is  due  to  an 
oversimplified  physical  picture  of  the  phenomena  or  to  too  crude 
mathematical  treatment  within  the  frame  of  the  physical  assumptions. 

For  this  reason  we  attacked  certain  key  problems  with  full 
mathematical  rigor.   This,  of  course,  involved  a  great  deal  of  work 
which  at  first  may  have  seemed  unwarranted  in  view  of  the  practical 
objectives  of  this  project*   Fowever,  our  efforts  In  this  direction 
paid  off  handsomely.   By  obtaining  a  rigorous  mathematical  result, 
we  are  able  to  judge  the  validity  of  the  physical  model.   Moreover, 
these  considerations  led  to  the  discovery  of  powerful  approximation 
methods  which  can  be  applied  to  the  problems  treated  here  and  to 
many  other  problems  of  a  similar  nature. 
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